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Part  7:  Some  Electrormgnetic  Problems  Relevant  to  the  EEG 
Paradigm 

Motivation  and  Background 

It  is  well  known  that  living  tissues  generate  electrical  activities.  The  electrical  signals 
produced  by  the  heart,  known  as  the  electrocardiogram,  the  ones  generated  by  muscles,  called 
electromyogram,  and  the  one  by  the  retina,  named  electroretinogram,  have  been  studied  by 
physicians,  neuroscientists1,2,3,  and  medical  researchers.  The  brain  also  produces  electrical 
activity,  and  the  study  of  such  activity,  i.e.,  electroencephalography  (EEG),  has  been  the  subject 
of  research  for  neuroscientists,  medical  researchers  and  clinicians.4,5  Since  these  electrical 
signals  are  due  to  the  stimulus-induced  potentials  and  currents  in  the  brain  tissues,  understanding 
and  analyzing  the  electromagnetic  characteristics  of  the  potential  and  currents  in  inhomogeneous 
dissipative  media  becomes  of  paramount  importance.  We  have  studied  and  analyzed  some 
selected  electromagnetic  problems  that  relate  to  the  EEG  paradigm. 


(1)  Electromagnetic  model  for  the  interpretation  of  the  potential  distribution 
over  the  scalp  due  to  electrical  activity  inside  the  brain:  Quasi-Static  Approach 

From  the  electromagnetic  point  of  view  we  can  think  of  the  brain  as  an  inhomogeneous,  lossy 
dielectric  body,  with  complex  geometry.  The  electrical  activity  of  the  brain  is  represented  by  the 
currents  flowing  inside  it  that  is  caused  by  the  neuronal  activities.  We  are  interested  in  the 
electromagnetic  modeling  and  analysis  of  the  effects  of  those  currents  on  the  potential  that  can 
be  measured  on  the  scalp  by  EEG  sensors.  In  order  to  have  a  qualitative  understanding  it’s 
possible  to  use  simplified  models  either  for  the  whole  brain  or  of  some  specific  part  of  it.  A 
simplified  model  of  the  brain,  effective  as  far  as  the  potential  distribution  over  the  scalp  is 
concerned,  is  well  known  multi-shelled  (i.e.,  multilayered  dielectric)  sphere  as  depicted  in  the 
Figure  1 .  In  this  way,  the  shape  of  the  brain  is  approximated  by  a  three-shell  spherical  volume 
with  the  different  complex  permittivities  in  the  various  layers4. 


1  A.  L.  Lyubarsky,  and  E.  N.  Pugh,  Jr.,  “The  recovery  phase  of  the  murine  rod  photoresponse  reconstructed  from 
electroretinographic  recordings”  Journal  of  Neuroscience,  vol.  16,  No.  2,  pp.  563-571,  1996. 

2  A.  L.  Lyubarsky,  B.  Falsini,  M.  E.  Pennesi,  P.  Valentini,  and  E.  N.  Pugh,  Jr.,  “UV-  and  midwave-sensitive  cone-driven  retinal 
responses  of  the  mouse:  a  possible  phenotype  for  co-expression  of  cone  photopigments”,  Journal  of_ Neuroscience  vol.  19,  No.  1, 
pp.  442-455, 1999. 

3  E.  N.  Pugh,  Jr.,  B.  Falsini,  and  A.  L.  Lyubarsky,  (1998)  “The  origin  of  the  major  rod-  and  cone-driven  components  of  the 
rodent  electroretinogram,  and  the  effect  of  age  and  light-rearing  history  on  the  magnitudes  of  these  components”.  In  T.P. 
Williams  and  A.B.  Thistle  (Eds.),  Photostasis  and  Related  Topics,  pp.  93-128.  New  York:  Plenum,  1998. 

4  Paul.  L.  Nunez,  Electric  Fields  of, the  Brain:  The  Neurophvsics  oj_ EEG,  Oxford  University  Press,  1981. 

5  David  Regan,  Human  Brain  Electronh  vsioloev:  Evoked  Potentials  and  Evoked_  Magnetic  Fields  in  Science  and  Medicine. 
Elsevier,  New  York,  1989. 
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Figure  1 :  Three-shell  model  of  the  brain 

The  structure  depicted  in  Figure  1  may  be  analyzed  either  from  the  electrostatic  or  the 
electrodynamic  point  of  view.  Since  the  frequencies  involved  are  very  low,  we  can  first  consider 
the  quasi-statis  approach.  . 

In  this  model  the  three  regions  represent  brain,  skull  and  scalp  and  the  whole  structure  is 
surrounded  by  unbounded  free  space.  Each  region  has  a  different  conductivity,  and  the 
surrounding  medium  has  zero  conductivity.  We  use  this  model  to  estimate  the  potential  in  each 
of  the  three  regions  and  outside,  due  to  the  current  sources  located  at  a  point  inside  the  brain. 
We  want  to  estimate  the  potential  distribution  in  the  structure,  related  to  different  configurations 
of  the  sources.  The  electrostatic  problem  to  solve  is  the  one  of  a  current  source  inside  a 
dissipative  medium.  From  the  continuity  of  charge,  we  have 

1 


In  the  quasi-static  case,  (l)is  equal  to: 


V-JDO  (2) 

Furthermore,  we  assume  that  the  Ohm  law  is  applicable  in  our  medium,  so  we  can  express  the 
current  density  as  the  sum  of  an  impressed  current  and  a  homic  current: 

J  =  J^+<rE  (3) 

where  a  is  the  conductivity  of  the  medium.  For  a  “Point  Source”  we  mean  a  point  from  which 
an  impressed  current  flows  into  our  medium.  We  will  call  a  “negative  source”  a  current  sink.  If 
we  substitute  (3)  in  the  equation  (2)  we  obtain: 

-V  •  crE  =  V  •  Jimp  (4) 

In  the  quasi-static  approximation,  the  electric  field  may  be  expressed  as  the  negative  gradient  of 
a  scalar  potential: 


E  =  -VT 


(5) 
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So  we  have: 


*  i 


?’0W  =  V.|„  (6) 

Equation  (6)  has  the  same  form  of  the  Poisson  equation: 

V-fiW  =  p.  (7) 

Equation  (7)  is  identical  to  equation  (6)  if  we  replace  £  by  cr  and  pbyV  -J.  ,  so  as  pointed  out 

by  Ref  [4],  we  have  only  to  solve  a  standard  electrostatic  problem  and  then  make  the  substitution. 
For  the  case  of  a  point  source  of  current  in  a  homogeneous  medium  with  conductivity  <x  we 
have: 


'?('•)  =  (8) 

Anu  r 

where  I  is  the  total  current  flowing  from  the  source.  Often  a  current  dipole  is  defined  as  a 
source  and  a  sink  of  current  very  close  to  each  other  and  with  the  same  magnitude. 


Radial  and  Non-Radial  Current  Dipoles 

The  first  basic  problem  to  solve  is  that  of  a  radial  current  dipole.  From  the  electromagnetic  point 
of  view  we  must  solve  the  Poisson’s  Equation  in  an  inhomogeneous  medium.  (We  start  with  a 
current  dipole  rather  than  a  simple  current  source,  because  we  assumed  that  the  surrounding 
medium  has  a  zero  conductivity  and  thus  the  current  flow  has  to  “close”  to  a  sink  in  a  region  of 
interest.  In  other  words  in  our  quasi-static  analysis  we  have  to  put  a  source  and  a  sink  of  the 
same  magnitude  in  order  to  “close  the  circuit”,  and  avoid  a  (time  dependent)  accumulation  of 
electric  charge  at  the  interface  between  conducting  and  non-conducting  (outside)  media.) 
Mathematically,  we  can  still  find  a  solution  for  a  source  and  one  for  a  sink,  and  then  superimpose 
them  before  we  impose  that  the  conductivity  of  the  surrounding  medium  is  equal  to  zero.  Let’s 
now  define  the  geometry  of  the  problem: 


Figure  2.  A  radial  current  dipole  in  the  three-shell  model 
We  will  use  the  following  names  for  the  different  regions: 

•  Region  1  (Brain):  r  <  Rx ,  conductivity  cr, 

•  Region  2  (Skull):  <  r  <  R2 ,  conductivity  cr2 

•  Region  3  (Scalp):  R2<r<R3,  conductivity  <r3 

•  Region  4  (Free  Space):  r>R3,  conductivity  cr4  =  0 


For  the  current  dipole  we  assume  the  following  parameters: 

•  Source  Position  r  =  rd  +  ^,  Source  Magnitude  I 

•  Sink  Position  r  =  rd-^.  Sink  Magnitude  I 

We  need  to  solve  the  Poisson’s  Equation  in  Region  1  with  the  source  and  with  the  sink,  and  the 
Laplace’s  equation  in  the  other  regions.  We  can  express  the  potential  in  Regionl  as  the  potential 
due  to  the  current  Source  in  free  space  plus  potential  due  to  the  “scattered”  field  at  the  interface: 


'F1='FS+VF 


si 


(9) 


where  we  put: 
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^=- 


Because  of  the  azimuthal  symmetry  of  the  problem  the  other  term,  which  is  a  solution  of  the 
homogeneous  equation,  is: 

(11) 

«=0 

Where  Pn  (cos  0)  is  the  Legendre  Polynomial  of  degree  n .  In  order  the  potential  to  be  finite  it 
should  be: 


K=  o 


In  Region  2  we  have: 


'Vs>=t1(v’+b1,r-')P„(cose) 


In  Region  3  we  have: 


=  E  («*»>•"  +b3„r~n~1)Pn(cos0 ) 


and  finally  in  Region  4: 


=  ZK  nr"  +b4nr-"-')Pn(c  os0) 


But  in  order  the  potential  to  vanish  at  infinity  it  must  be: 


a4„  =0 


In  order  to  impose  the  boundary  conditions,  we  can  expand  the  potential  (10)  due  to  the  Current 
Source  using  the  orthogonality  properties  of  the  Legendre  Polynomials,  as  follows: 


vp  =. 
T  s 


r-  r.+ 


m  mm  r,K+- 

^T=7 £7X7  r  , - ^P"(cos0)  (17) 

d)  Ancrx  „=0  f  f  (  d\T) 


max  r,\  rd+- 
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Now  we  impose  the  boundary  conditions  at  each  interface.  At  the  first  interface  we  have: 


cWs(r)  9¥sl(r) 

dr  n  dr 


WM(r) 


At  the  second  interface: 


*«(*»)  = 

.  ^wi  __  w«wi 


And  at  the  third  one: 


_  a^Ml  .  w-.wl 


Based  on  what  we  said  before  we  cannot  set  cr4  =  0  yet.  (One  can  show  that  if  we  set  cr4  =  0  at 
this  stage,  the  linear  system  has  an  indeterminate  solution.  This  can  be  explained  as  follows: 
Let’s  assume  that  cr4  =  0 ,  in  this  case  the  second  equation  in  (20)becomes: 


So  all  coefficients  other  than  a2  are  equal  to  zero,  so: 


If  we  put  (22)  in  (19)  we  obtain 


*«(*,)■ =«3 


And  finally  if  we  put  (23)  in  (18)  we  obtain: 
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(24) 


'Ps(i?I)  +  ^1(i?1)  =  «3 


ws(r) 


dr 


+ 


3^,(0 


dr 


=  0 


And  from  (24)  we  have 


(25) 


and  we  cannot  get  a  unique  solution.) 

We  can  solve  the  same  problem  for  the  negative  source  and  add  the  result  to  the  one  obtained  for 
the  positive  source.  Only  at  this  point  we  can  set  cr4  =  0  to  get  the  final  result  for  the  potential 

distribution  due  to  the  current  dipole.  In  Figure  3  we  show  the  distribution  of  electric  potential 
due  to  a  radial  current  dipole.  In  the  figure  we  can  see  the  equipotential  lines.  We  used  the 
following  numerical  parameters: 

Region  1  (Brain):  radius  =  7.65  cm  cr,^— — 

1  350 

Region  2  (Skull):  thickness  =  0.8  cm  cr  = — - —  (26) 

2  20000 

Region  3  (Scalp):  thickness  =  0.2  cm 
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Figure  3.  Potential  distributions  due  a  current  dipole  in  three-shell  model 
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Using  the  same  technique  we  can  analyze  the  case  depicted  in  Figure  5  of  non  radial  current 
dipole. 


Figure  4:  A  non-radial  current  dipole  in  three-shell  model 

In  this  case,  we  no  longer  have  the  azimuthal  symmetry  of  the  previous  case,  but  we  can  still 
separately  calculate  the  contribution  of  the  positive  and  the  negative  source  and  then  we 
superimpose  the  effects.  .  The  result  of  this  simulation  is  shown  in  Figure  5.  It  is  worth 
highlighting  from  the  Figures  3  and  5  that  a  distribution  of  potential  could  be  measured  even  at  a 
point  away  (i.e.,  non-contact)  from  the  scalp  surface. 
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Modeling  of  a  Probe:  Perturbed  Solution 


How  is  the  potential  distribution  perturbed  by  the  presence  of  a  conductive  small  probe  inside  the 
structure  or  on  the  surface?  The  configuration  that  we  have  analyze  is  shown  in  Figure  6. 


Figure  6 


We  assume  that  the  radius  of  the  spherical  probe  is  rp  and  its  conductivity  iscrp.  We  use  the 

perturbation  approach.  First,  we  calculate  the  potential  on  the  structure  without  the  probe  and 
evaluate  the  electric  field  at  the  point  where  the  probe  is  to  be  located.  Since  the  probe  is  very 
small  we  assume  that  the  incident  field  is  “uniform”  over  the  region  that  is  to  be  occupied  by  the 
probe,  and  is  equal  to  the  unperturbed  field  in  the  center  of  the  probe.  Let’s  denote  the 
unperturbed  potential  as  vFo(r,0,^) .  The  electric  field  will  be: 

E  fr,G,(p)  =  -'Wfr,0,(p)  (27) 

If  the  position  of  the  probe  is{ro,0o,<po} ,  the  electric  field  that  excites  the  probe  is: 

E0  =Eo(ro,0o,po)  (28) 


Due  to  these  assumptions,  this  field  induces  a  uniform  field  inside  the  probe,  equal  to: 
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(33) 


p  a  n+2cr, 


En 


rp  (r-r0)-E0 


r-rn 


r~r0  E0 


This  is  a  first-order  approximate  of  the  perturbation  due  to  the  presence  of  the  probe.  This  can 
be  further  improved  as  follows:  After  the  introduction  of  the  perturbed  field  in  (33)  due  to  the 
probe,  the  boundary  conditions  are  no  more  satisfied,  so  we  have  to  solve  again  the  problem  of 
Poisson  equation,  with  the  induced  dipole  on  the  probe  as  a  secondary  source.  We  will  then 
superimpose  the  perturbed  potential  with  the  unperturbed  potential.  We  assume  that  the  induced 
dipole  on  the  probe  is  due  to  a  positive  and  a  negative  source  of  current  close  to  each  other  in 
such  a  way  that  the  induced  current  dipole  is  given  as: 


P  =  4^i  —  P  rpEo  =  4»oi/pL 
cr  +2C7, 


(34) 


which  is: 


<7  -  <7,  , 

-J*—Lr%=Ip  L 


cr  +20-, 


(35) 


Our  constraint  is  on  the  dipole  moment,  so  we  have  a  degree  of  freedom  in  choosing  I p  and|L| 
We  arbitrarily  choose  |L|  to  be  equal  to  the  diameter  of  the  probe  2 rp ,  and  then  we  have: 


I  r± 
p  a p  +  2  cr,  2 

L  =  2r 

'  Ej 


lEo 


(36) 


In  this  way  the  two  sources  are  given  by: 


(37) 


V(r)  =  // 


r- 


ro+rP 


En 


E, 


oU 


Ur)-W 


r- 


«iy 


As  in  the  unperturbed  potential,  we  can  now  solve  two  separate  problems  with  azimuthal 
symmetry  and  then  superimpose  the  solutions  to  have  the  total  perturbation  potential.  In  the 
Figures  8  and  9,  the  green  “dot”  shows  the  position  of  the  probe. 

As  we  can  see,  there  is  almost  no  change  in  the  potential  provided  the  probe  is  small. 


(2)  Electromagnetic  model  for  the  interpretation  of  the  potential  distribution 
over  the  scalp  due  to  electrical  activity  inside  the  brain:  Electrodynamic 
Approach 

We  can  also  analyze  this  structure  using  the  electrodynamical  approach.  In  this  case  we 
substitute  the  current  dipole  used  in  the  quasi-static  analysis  by  an  electromagnetic  (fipole.  From 
the  electromagnetic  view  point,  we  have  to  solve  the  problem  of  the  radiation  by  an  electric 
dipole  in  a  multi-layered  medium  with  different  finite  conductivities. 

Tai  and  Collin  in  their  2000  paper  in  the  IEEE  Transactions  on  Antennas  and  Propagation  (Vol. 
48,  No.  10,  October  2000,  pp.  1501-1506)6  discussed  an  interesting  problem  dealing  with 
radiation  of  a  Hertzian  dipole  in  a  conducting  medium.  They  concluded  that  if  a  dipole  is  totally 
immersed  in  a  finitely  conductive  medium,  for  a  given  dipole  moment  the  total  radiating  power 
can  theoretically  be  infinite.  They  mention  that  this  is  due  to  the  dissipation  of  energy  in  the 
immediate  vicinity  of  such  a  dipole  where  the  field  can  in  principle  be  infinitely  strong.  Then 
they  analyzed  the  case  of  an  insulated  dipole  in  a  conducting  medium,  in  which  the  dipole  is 
insulated  by  a  lossless  dielectric  sphere  and  then  it  is  immersed  in  a  finitely  conducting  medium 
of  infinite  extent.  Their  results  show  that  in  this  case  the  radiated  power  is  finite.  This  problem 
has  direct  relevance  to  the  EEG  problems  in  which  induced  dipoles  in  the  brain  tissues  are 
effectively  dipoles  in  a  conducting  medium,  resulting  in  the  volume  current  in  this  medium. 
Specifically,  it  is  known  that  the  pyramidal  cells,  the  elongated  primary  processing  neurons  of 
the  visual  cortex,  necessarily  act  effectively  as  “current  dipoles”  because  of  the  spatial 
distribution  of  their  membrane  channels,  so  that  their  local  “sources”  and  “sinks”  lie  at  some 
distance  from  one  another  along  a  line  perpendicular  to  the  cortex.  Therefore,  understanding  the 
behavior  of  such  dipoles  in  finitely  conducting  medium  is  important  in  the  EEG  study.  Here  we 
have  explored  the  analysis  and  modeling  of  fields  due  to  insulated  and  small  dipoles  in  a 
conducting  sphere  of  layers  with  different  conductivities:  We  refer  to  the  geometry  in  Figure  8 
where  in  general  the  insulated  dipole  is  off-centered. 


®  C.  T.  Tai  and  R.  E.  Collin,  “Radiation  of  a  Hertzian  Dipole  Immersed  in  a  Dissipative  Medium,”  IEEE  Tmnsactigns  on 
Antermas  SSi Propagation,  Vol.  48,  No.  10,  October  2000,  pp.  1501-1506 
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Figure  8.  Radiation  of  a  dipole  in  three-shell  model 


To  solve  this  problem  we  expand  the  electric  and  the  magnetic  fields  in  terms  of  Spherical 
Vector  Wave  Harmonics7  in  each  region  and  with  respect  to  the  most  convenient  reference 
system  between {r,0,(p) and {r',0',(p'} .  Spherical  Vector  Wave  Harmonics  are  defined  as  follows: 

M»,«  (r  ,0,  <p)  =  V  unm  (r  ,0,  (p)xr 

N„(r,«,p)  =  |vxM„(r,S,p)  (38) 

v.,(r,e,V)-j „(kr)p;(cosey” 

N<“i(r,«,f>)  =  ivxMS(r,0,«,)  (39) 

k 

In  region  1  we  expand  the  fields  with  respect  to  {r,0,<p} : 

E ,(r,0, •?)  =  ££  (40) 

77=1  m=-n 

Ufr,0,<p)  =  ^^J]Pl(n>m)Mnm(r,0,<p)-t{)l(n,m)Nnm(r,e,(p)  (41) 

7l  n=l  m=-n 

1  A.  Stratton,  u Electromagnetic  Theory”,  McGraw-Hill,  1941 
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Due  to  the  presence  of  the  electric  dipole  in  the  region  1  we  have  to  write  the  fields  in  the 
following  way: 


E,  (r,0,(p)  =  (1,  0)N{Jq  (r  ,0,<p)  +  /}(2)  (1, 0)N|20)  (r  ,9,<p)  + 

+P\  (1»  -1)N,  (r,  0,<p)  +  Pl  (1, 1)N,  (r,  0,  <p)  + 

+Z  Z  ^(«>W)N^(^^)  +  Z  Y  e>,m)M,,m(r,0,p) 


n= 2  m=-n 


m=-n 


(42) 


H,  (r,0,<p)  =  ^[P^\\,0)M%{r,e,(p)  +  O)Mg(r,0,p)]H 

+ [i>  (1,  -  OM,.,  (r,0,p) + Px  (U)M,  (r,0,<p)]  + 

+-Z  Z  + £  X  Q(n,m)NB>,0,p) 

*7l  w-2  m=-/t  «=1  m=-n 


(43) 


In  region  2  we  expand  the  fields  with  respect  to  {r,0,<p} : 


E2(r,0,<P)  =  f;  £  (r,tf,p)+Qf  (r, «,«.)  + 

w=l  m=-« 


«=1  m-~n 


H2(i-,0,*>)-— X  ^  P2<‘»(«,m)M™  (r,0,p)-t£<Vm)N®  (r, £,,»)  + 

^7 2  fl=l  m=-rt 

+— Z  Z  (/-,£, (/•,»,?>) 

Vl  n=l  m=-« 


(44) 


(45) 


In  region  2  we  also  expand  the  fields  with  respect  to  {r\ 6',  cp'} : 


E 2(r\ 0  >') =£X4V  «)N®  (/•',  ff,  <p')  +P®  („,  /»)M®  (/■',<?>') + 


/i— 1  m=~n 


+Z  Z  42)(«»ra)NS(^'^'>^')+^22)(«»w)M(^(r',^>') 


«=1  m=-n 


(46) 


^7 2  «=1  m=-n 

+-Lfl  Z  4,)(».'»)mS(>-,,p>')+s?)(».'»)nS(/,  «>') 

02 


(47) 
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In  region  3  and  in  region  4  we  expand  the  fields  with  respect  to  {r',0',cp'} : 


E 3(r',0>’)  =  £  £  4,,(n,m)N®  (I'.fl'.rt  +4"(n,m)M®  (r',0» 

«=1  m=-n 


(48) 


n= 1 


£  4"(«,»>)M®  (/•',£>', p')+4”(«.»<)N®  (r’AfO- 

Vi 

+— £  £  4!,(n, «>')+«», m)N®  (/,<?>') 

^3  /i=l  m=-« 


(49) 


In  region  4  we  impose  also  the  radiation  condition: 


E «(/, «>')  =  £  £  42>(n,m)N®  (r',e',/)+42»(n,m)M®  (r’,0>')  (50) 


«=1  m--n 


H4(r',#',  <p')  =  —  £  £  42>(»i,m)MS(/,0',«>')+S<2,(n,m)NS(r',e',«)O  (51) 

^?4  «=1  m——n 


By  imposing  the  boundary  conditions,  we  can  find  the  unknown  coefficients.  Obviously,  we  can 
easily  impose  the  boundary  conditions  for  fields  expressed  in  the  same  reference  system,  so  we 
can  impose  the  boundary  conditions  at  the  interface  between  regions  1  and  2  using  the  reference 
system  {r,  0,  <p) .  At  the  interface  between  regions  2  and  3,  we  use  the  reference  system  {r',  & ,  <p'} 
Then  we  have  to  impose  the  Source  Condition: 


^(I)  (1,0)  +  ^(2)  (1,0)  =  -^- 

4  n 


(52) 


Finally  we  have  to  relate  the  field  in  region  2  expressed  in  both  reference  systems  by  using  the 
addition  theorem  for  spherical  vector  wave  functions. 

Let’s  consider  the  electric  field  in  the  region  2  written  in  both  reference  systems: 

E2(r,0,?>)  =  £  £  lf(^»I)N2,(r,«,rtlifl(i>,»)M®M,f)t 

n= 1  m=-n 

+£  £  P?\nMKZ(rA9)^?\n,m)M™(r,e,<p) 

n=l  m=-n 

E,  (r\  ff,  <p')  =  £  £  4‘>  (n,  m)N®  (r’.  +B f  („,  m)M®  (r\ 0>')  + 

/7=1  m~-n 


(53) 


n- 1  m=-n 
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The  Addition  Theorem  for  the  spherical  vector  wave  functions8  states: 


V~\  fl-~V 

V=\  fl=-v 


(54) 


where: 


X^\n,m,v,M)  =  (-irY^,rr\-M,y\p,P)c(n,v,p)h{;\kr0)P;-^^ 

P 

ya)  (n,  /n,  v,  p)  =  (-l)/,+1  a(m,  n|  -p,v\p,p-  l)d(n,v,  p)hf  (kr0)P”~,i  (cos  0o)ei(m~f,)<p,> 


c(n,v,p)  =  i 


_  . y+p-n  [2v(v  +  l)(2v+l)  +  (v +!)(«  —  v  +  p +l)(n  +  v - p)—v(v -n  +  p  +  V)(n  +  v  +  p  +  Tj\ 


d(n,v,  p)  =  i 


2v(v+l) 

.y+P-„  (2v  +  l)yl(n  +  v  +  p  +  l)(y-n  +  p)(n-v  +  p)(n  +  v-p  +  Vj 


2v(v  + 1) 


(55) 


Using  that  theorem  we  can  write  the  field  in  region  2  as: 
E2  (r  E2  (/,  O',  cp')  = 


oo  n 


=E  E  3V®)  E  E  *(V™,v,//)n^ 

^  v=l  //=-v 


11=1  m=-n 
oo  n 


+ 


+E  E  CgV™)  E  E  (r',0', <f,')+X"\n,m,v,tW'-»y,e\<p") 

n= 1  m--n  v~\  p~~v  j 

+E  E  P^mit  X  A(2)(«,m,v,//)N(^^^  1  + 

n=l  m=-n  y  i/=l  p-~v  J 

+E  E  efw/E  E  *(2)(«,™,v^)m^ 

«=1  m~~n  y  i/=l  fi=-v 

(56) 

At  this  point  we  change  the  order  of  the  summations  in  (56)  and  we  obtain: 


8  O.R.  Cruzan,  “ Translational  Addition  Theorems  for  Spherical  Vector  Wave  Functions”,  Q.  Appl.  Math.,  vol.  20,  pp  33-40, 
1962 
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*  i  * 


E2(r',0\<p')  = 


=EE 


2  2  ^20) («> w)^0) 0, w, v, /i)  +  (n, m)Ym (n,m,v,/j) 


v=\  n=-v[_  «=1  n=-n 
co  v  P  co  n 


v/=l  fj=-v  L  «=1  m=-« 


uj  K 

EE 


2  is  P?\n,m)X(2\n,m,v,v)  +  Q2\n,m)Y('2\n,m,v,iu) 


v- 1  //=-vL«=l  n=-rt 


M®(r',0>')  + 
N®(/,0>')  + 


W  K 

+EE 


J  2  -P2(2)(M’W)I"(2)(«»W»l/^)  +  022)(^W)^(2)(^WJV/>//) 


k=1  fJ=-v  L  w=l  m--n 


M"i  (r’.e'.fO 


(57) 


At  this  point  we  can  compare  (53)  with  (57)  and  we  obtain  the  following  relations  between 
coefficients: 


*■ 

41) (v>/i)  =  EE  Pl1\n,m)Xm(n,m,v,ju)  +  Q2\n,m)Ym(n,m,v,/j;) 

n- 1  n=-/i 

#2 ,)0'»/0  =  2]  £  P2(1)(n,m)F(1)(n,w,v,//)  +  ^1)(n,m)A(1)(n,/w,v,//) 

"I’T  (58> 

42)  O'* M)  =  a  P22) m)^(2) («, v, //)  +  022) («, w)F(2) («,  m,  v,  //) 

n=l  n=-n 

Bi\v’P)  =  YjYj  Pl2\n,m)Y(2\n,m,v,{i)+Q<?\n,m)Xi2\n,m,v,{i) 

_  /I=l  /»=-« 

Finally  using  the  conditions  (58)  we  can  find  all  of  the  unknown  coefficients. 

It  is  important  to  mention  that  when  we  impose  the  boundary  conditions  at  the  interfaces  1-2  and 
2-3,  because  of  the  orthogonality  properties  of  the  spherical  vector  wave  functions,  we  see  that 
all  of  the  equations  are  decoupled,  so  we  can  find  immediately  a  solution.  The  only  equations  we 
have  to  solve  are  the  conditions  (58).  In  other  words  the  real  number  of  unknowns  depends  only 
on  the  regions  with  non  concentric  boundaries,  in  our  case  the  region  2. 

The  geometry  becomes  much  simplified  if  the  insulated  dipole  is  located  at  center  of  the  sphere9, 
as  shown  in  Figure  9 


9  C.T.  Tai,  R.E.  Gollin,  “Radiation  of  a  Hertzian  Dipole  Immersed  in  a  Dissipative  Medium”,  IEEE  Trans.  Antennas  and 
Propagation,  vol.  48,  n.  10,  October  2000 
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Figure  9 


In  this  case  we  have  3  concentric  dielectric  spheres  with  an  electric  dipole  located  in  the  center. 

In  order  to  find  electric  and  magnetic  fields  everywhere  we  write  the  Helmholtz’s  equation  for 
the  magnetic  vector  potential  in  the  four  regions: 


v24z(r)  +  tfAX!(r)  =  -jjJS(r)  r  <rx  (59) 

V2A2z(r)  +  k22A2z(r)  =  0  rx<r<r2  (60) 

V24z  (r)  +  £324z  (r)  =  0  r2<r<r2  (61) 

V24  z(r)  +  k]A^{  r)  =  0  r>r3  (62) 


Due  to  symmetry,  the  solutions  are  written  as: 


A,.  = 


f  e~Jk'r  e 

a, - +  b, 


jk,r  \ 


'i 


V 


I  =  1,2, 3, 4 


(63) 


In  order  to  find  8  unknown  coefficients,  we  have  to  impose  boundary  conditions  at  each  interface 
(6  equations),  the  radiation  condition  (1  equation)  and  the  source  condition  (1  equation). 

The  radiation  condition  states  that  in  the  unbounded  region,  for  time-dependence  eJmt ,  we  get : 


*4=0 


(64) 
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Jl 


For  the  source  condition  we  consider  the  inhomogeneous  Helmholtz’s  equation  (59)  and  we 
integrate  both  sides  in  the  volume  V  of  a  sphere  of  radius  a  -»  0 : 

fvVA,,(r)dr+k!  J4,(r>fr  = -/a/  r<r,  (65) 

V  V 

The  second  integral  in  the  left  side  is  equal  to  zero  for  a  — »  0 .  We  apply  the  divergence  theorem: 


2  It  JT 

1) 

0  0 


dAlz(r  ,0,<p) 


dr 


a2  sin  0d0d<p  =  -fj,J 


r  <r. 


(66) 


Then  we  put  (63)  into  (66)  and  we  obtain  the  source  condition: 

a, +4,  =7^  (67) 

4  n 

For  the  other  boundary  conditions  we  have  to  calculate  electric  and  magnetic  fields  in  each 
region: 


-VxA, 


V(V-  A,) 


~M A, 


(68) 

(69) 


And  then  we  can  impose  the  boundary  conditions: 


rxEj 

rxHi 

r  xE2| 

r  xH2 

rxE3| 

fxH3 


=  r  xE2| 
=  rxH2| 
=  rxE,| 
=  r  xH3 


(70) 

(71) 

(72) 


We  have  to  solve  the  system  of  equations  given  by  (64),  (67),  (71)  and  (72).  Once  we  have  the 
solution  of  the  system  we  can  plot  the  behavior  of  the  fields.  In  figure  10,  we  show  the  behavior 
of  the  absolute  value  of  the  transverse  electric  field. 
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Figure  13.  Magnitude  of  the  Transverse  electric  field  of  the  insulated  dipole  in  finite-radius 

sphere 


In  Figure  11,  we  show  Poynting  Vector  (multiplied  byr2),  i.e.,  the  power  radiated  per  unit  solid 
angle,  in  order  to  see  the  effect  of  dissipation  due  to  the  loss: 


Figure  14 

In  the  outer  region,  without  loss,  the  function  is  constant,  as  expected. 


* 


(3)  Potential  of  a  Small  Dipole  in  a  Dissipative  Medium  with  Certain  Inhomogeneous 
Conductivity 

In  our  previous  study10,  we  found  that  the  sagittal  fissure  in  the  brain  can  act  as  a  “major  shunt” 
for  electrical  current  flow,  providing  evidence  that  the  commonly  used  3-sheel  model  of 
conductivity  may  not  be  adequate  in  certain  conditions.  This  is  due  to  the  fact  that  the 
conductivity  in  the  region  of  sagittal  fissure,  which  is  filled  with  the  cerebrospinal  fluid  (CSF),  is 
different  from  that  of  the  other  inner  parts  of  the  brain.  In  that  work,  it  was  shown  that  a  pair  of 
visually-evoked  induced  “dipole  currents”,  i.e.,  excited  neurons  in  the  visual  cortex,  oriented 
~180  degrees  opposite  and  located  across  the  sagittal  fissure  produce  the  scalp  potential 
distributions  nearly  identical  in  amplitude  and  phase,  whereas  the  3 -shell  model  suggested  a 
mirror-image  symmetry  across  the  midline  for  this  potential  distribution.  Furthermore,  the  scalp 
potential  around  the  region  of  projection  of  the  sagittal  fissure  onto  the  scalp  had  maximum 
amplitude,  which  is  again  different  from  what  the  3-shell  model  predicted.  Here  we  have 
explored  this  issue  theoretically  by  analyzing  certain  electromagnetic  model  of  a  current  dipole 
in  a  medium  with  certain  inhomogeneous  conductivity.  We  have  studied  theoretically  the 
potential  distribution  of  a  current  dipole  in  a  conducting  medium  with  uniform  conductivity 
everywhere,  except  in  a  wedge  (or  a  cone)  region  having  a  different  conductivity  (Fig.  12), 
essentially  modeling  sagittal  fissure  shunt.  The  current  dipole  can  be  located  at  some  arbitrary 
point.  We  utilized  a  quasi-static  model  for  this  problem  and  have  expressed  the  electric 
potentials  everywhere. 


10  E.  N.  Pugh,  Jr.,  J.  B.  Demb,  B.  Alterman,  C.  Daniele,  L.  Daniele,  H.  Baseler,  N.  Engheta,  and  B.  Wandell,  "Visual  potentials 
evoked  by  small  oscillating  sources  located  by  fMRI:  the  sagittal  fissure  is  a  major  shunt,"  a  poster  presented  in  the  29th  Annual 
Meeting  of_ the  Society,  (or  Neuroscience .  in  Miami  Beach,  Florida,  October  23-28, 1999. 
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